Quantum lump dynamics on the two-sphere 
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<JJ ■ Abstract 

It is well known that the low-energy classical dynamics of solitons of Bogomol'nyi 
type is well approximated by geodesic motion in M n , the moduli space of static tl- 
solitons. There is an obvious quantization of this dynamics wherein the wavefunction 
%j) : M n — > C evolves according to the Hamiltonian Hq = ^ A, where A is the Laplacian on 
' M n . Born-Oppenheimer reduction of analogous mechanical systems suggests, however, 

^sO ■ that this simple Hamiltonian should receive corrections including k, the scalar curvature 

Tj" ', of M n , and % ', the n-soliton Casimir energy, which are usually difficult or impossible to 

compute, and whose effect on the energy spectrum is unknown. This paper analyzes the 
spectra of Hq and two corrections to it suggested by work of Moss and Shiiki, namely 
H\ = Hq + -iK and H2 = H\ + in the simple but nontrivial case of a single CP 1 
lump moving on the two-sphere. Here Mi = Rati, a noncompact kahler 6-manifold 
invariant under an SO(3) x 50(3) action, whose geometry is well understood. The 
symmetry gives rise to two conserved angular momenta, spin and isospin. By exploiting 
the diffeomorphism Rati — TSO(3), a hidden isometry of Rati is found which implies 
that all three energy spectra are symmetric under spin-isospin interchange. The Casimir 
energy is found exactly on a SO(3) submanifold of Rati, using standard results from 
harmonic map theory and zeta function regularization, and approximated numerically 
on the rest of Rati. The lowest 19 eigenvalues of Hi are found, and their spin-isospin 
and parity compared for i = 0, 1, 2. It is found that the curvature corrections in Hi lead 
to a qualitatively unchanged low-level spectrum while the Casimir energy in H2 leads to 
significant changes. The scaling behaviour of the spectra under changes in the radii of 
the domain and target spheres is analyzed, and it is found that the disparity between 
the spectra of Hi and H2 is reduced when the target sphere is made smaller. 
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1 Introduction 



Many field theories arising naturally in theoretical high energy physics may be said to be 
of Bogomol'nyi type. For such theories there is a topological lower bound on the energy 
of field configurations, and this bound is attained only by solutions of a first order "self- 
duality" equation, the so-called solitons of the theory. The solitons are stable by virtue of 
their energy-minimizing property, and are generically spatially localized lumps of energy with 
strongly particle-like characteristics. Examples are magnetic monopoles, abelian Higgs vortices 
and sigma model lumps. There is a well-developed geometric framework for studying the 
classical low-energy dynamics of such solitons, proposed originally for monopoles by Manton 
[22], called the geodesic approximation. Here, n-soliton trajectories are approximated by 
geodesies in the moduli space M n of static n-solitons, with respect to the metric induced by the 
kinetic energy functional of the field theory (usually called the L 2 metric). This approach has 
been extremely fruitful, providing both important qualitative insight into topological soliton 
dynamics and good agreement with numerical analysis of the full field theory. In the case of 
vortices and monopoles, the geodesic approximation is backed by rigorous analysis [38, 39]. 
For a comprehensive review, see [23]. 

When one comes to quantize the low-energy dynamics of such solitons, two approaches 
are possible. Since geodesic motion on M n captures the classical soliton dynamics so well, it 
is natural simply to quantize that [13]. Then the quantum n-soliton state is specified by a 
wavefunction ip : M n — > C, evolving subject to the Hamiltonian 

H geo = ±A, (1.1) 

where A is the Hodge Laplacian on M n . This has the virtue of simplicity, but it ignores 
the degrees of freedom normal to the moduli space completely. An alternative is to make 
a low energy reduction of the full quantum field theory by means of the Born-Oppenheimer 
approximation. This has been carried out for sigma model lumps by Moss and Shiiki [28]. 
Arguing by analogy with finite dimensional mechanical systems, they find that, once again, 
the low energy quantum dynamics of n solitons can be described by a wavefunction on M n , 
but that the Hamiltonian is 

H BO = \^ + \n-\\\k\\ 2 + U + ■■■ , (1.2) 

where k is the scalar curvature of M n , k is the mean curvature of the embedding of M n into the 
(infinite dimensional) field configuration space and U is a potential on M n incorporating the 
residual effects of the normal modes. This potential is rather complicated, but its dominant 
term is the n-soliton Casimir energy, that is, the total zero-point energy of the normal modes 
to the static n-soliton. 

The aim of the current paper is to compare the spectra of H geo and Hbo, to determine to 
what extent the extra terms in (1.2) change the quantum energy spectrum. Of course, one 
expects the numerical values of the energy eigenvalues to change to some extent. In itself, 
this is not particularly important. We will be more interested in the qualitative features of 
the energy spectrum. For example, on a noncompact moduli space, one could easily imagine 

that H geo 

may have no bound states, while Hbo does. In the extreme case, one could find 
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that the spectrum of H geo is continuous while that of Hbo is discrete. Clearly, H geo would be 
a disastrously bad approximation in this case. Other, more refined, qualitative features can 
also be compared, for example, the dimension and symmetry properties of the eigenspaces of 
H geo and Hbo, ordered by energy. 

We address this issue for a specific example, namely a single CP 1 lump moving on the 
two sphere (spacetime S 2 x K). This choice has several mathematical advantages, making it 
simple enough to be tractable, but not too simple to provide a nontrivial test of the truncation 
H geo . First, the static n-solitons are explicitly known - they are rational maps, so M n = Rat n . 
Second, by choosing physical space to be compact (S 2 here) rather than M 2 (as for the model 
on Minkowski space M 2+1 ), we ensure that the L 2 metric on Rat n is well defined (there are 
no non-normalizable zero modes). Third, we may think of the L 2 metric (formally) as the 
induced metric on Rat n as a complex submanifold of (CP 1 ) 5 , the infinite dimensional space 
of maps S 2 — > CP 1 . This latter is (formally) an infinite product of kahler manifolds, so is 
kahler, so, as noted by Moss and Shiiki [28], the very awkward extrinsic curvature term is 
absent, ||fc|| 2 = 0. Fourth, the L 2 metric on Rati is highly symmetric, so the L 2 metric and its 
scalar curvature can be computed explicitly. On the other hand, Rati is not too symmetric: 
the scalar curvature and Casimir energy are non-constant functions of a single variable, the 
lump width, so that (in contrast to the case of monopoles and vortices) even the n — 1 sector 
provides a nontrivial test. 

To be precise, we will compute the spectra of three approximations to Hbo, 

H = H geo = ±A, Hi = |A + |ac, h 2 = \A + \k + ^ } (1.3) 

where ^ is the Casimir energy. We will find that the spectrum of H can be thought of as 
a perturbed version of the spectrum of the Laplacian on CP 3 (equipped with the Fubini- 
Study metric). Quantum lumps possess two integer conserved angular momentum quantum 
numbers, which we call isospin k and spin t, associated with the rotational symmetries of 
the target space CP 1 and the domain S 2 respectively. We will show that for each choice of 
k,t, the spectral problem for Hi, i = 0, 1,2, reduces to a matrix Sturm-Liouville problem of 
dimension 2min{/c,t} + 1. We find that this problem is symmetric under interchange of k 
and t, so the energy spectra of H ,Hi,H 2 all possess this symmetry. A careful analysis of 
the boundary conditions for the Sturm-Liouville problem shows that the boundary conditions 
appropriate for Hi are the same as for H 2 , but different from H : including scalar curvature 
changes the boundary conditions. Nonetheless, the spectrum is discrete in all three cases. 
We have computed (numerically) the first 19 energy eigenvalues for each Hamiltonian and 
tabulated {k,t} for the eigenstates in order of increasing energy. The spectra of Hq and Hi 
are remarkably similar. The order of states is the same apart from two transpositions. The 
Casimir energy leads to significant changes in the spectrum. However, the effect of the Casimir 
energy becomes less important if the radius of the target space is decreased. 

2 The one-lump moduli space 

Identifying CP 1 = S 2 and using complex stereographic coordinates z = x + iy and W on 
domain and codomain respectively, we may identify a map : S 2 — > CP 1 = S 2 with a 
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complex function W{z). The kinetic and potential energy functionals of the CP 1 model on 
S 2 are then 



T = 1 -[ w ^ifft[iH) d ^ (21) 

2 7 (1 + |^|2)2 (1 + 1^)2)2' V 2 J (1 + | W |2)2 ^ 

respectively. Note that, in contrast with earlier work [3, 25, 36, 37] we have given both domain 
and codomain spheres the round metric of radius 1. This will allow us to make direct use of 
results in the harmonic maps literature when we come to compute the Casimir energy of a 
lump. We will consider how our results change when the radii of the domain and target spheres 
are altered in section 5.5. 

A theorem of Lichnerowicz [20] shows that, if <fi has topological degree n (assumed, without 
loss of generality, to be non- negative), then V > 47m, with equality if and only if is holomor- 
phic. So degree n holomorphic maps : S 2 — > S 2 minimize V within their homotopy class: 
these are the solitons of the model, and the "self-duality" equation is the Cauchy-Riemann 
equation. In terms of W and z, degree n holomorphic maps S 2 — > S 2 are rational maps of 
(algebraic) degree n, that is, 

W = aiZH n + ' ' ' + an+1 , (2.2) 

where aj are 2n + 2 complex constants, a\ and a n+ 2 do not both vanish, and the numerator and 
denominator have no common roots. One interprets this physically as a static superposition 
of n lumps. So the moduli space of static n- lumps is Rat n , the space of degree n rational maps. 

The L 2 metric on Rat n is defined by restricting the kinetic energy functional T to fields 
W(t, z) which at each fixed t are degree n rational maps. Explicitly, in the chart where a± ^ 0, 
one can define local complex coordinates qi = a^/ai, i = 2, . . . , In + 2 on Rat„. Then allowing 
<?2, • • • , <?2n+2 to vary with time, one substitutes 

Wit Z) = Z n + ^ n ~ 1 + --- + «r l+ l(t) 

qn+2(t)z n + q n+3 (t)z n ~ 1 + --- + q 2n+2 (t) 

into T, to obtain 

T _1 V _ f 4 OWdW Adxdy 

~ 2 2_^Hm^ ^ ~ J (1 + \W\ 2 ) 2 dq t dq 3 (1 + \z\ 2 y K } 

The L 2 metric is 7 = 2~2i j lijdqidqj. See [37] for a coordinate free definition of 7. The metric 7 
is manifestly Hermitian. In fact, it is kahler [32, 37]. It is known to be geodesically incomplete 
[33], so the classical geodesic approximation predicts lumps may collapse to infinitely narrow 
spikes in finite time. Numerical simulation and rigorous analysis (albeit on domain C) confirm 
that lump collapse can occur, though the geodesic approximation gets the fine detail of the 
singularity formation process wrong [7, 21, 31]. 

There is an isometric action of G = SO(3) x 5*0(3) on Rat n , induced by the natural SO(3) 
actions on the domain and target spheres. On Rat x this action has cohomogeneity 1 (generic 
G orbits have codimension 1), and, in fact, almost completely determines 7. Consequently, an 
explicit formula for 7 is known in this case, and the geometry is particularly well understood. 



4 



For n — 1, the no common roots condition on the rational map W(z) = (a\z + a2)/{a 3 z + a 4 ) 
is aia 4 — a 2 a 3 7^ 0, so we may identify each map with a projective equivalence class [L] of 
GL(2, C) matrices, 

(i\Z + a 2 



1 M . (2.5) 
a 3 z + a 4 |_ V ° 3 ° 4 / 

Hence Rati = PL(2,C). The action of G on Rati corresponds, under this identification, with 
the natural action of PU{2) x PU{2) on PL(2, C): 

([C/!],^]):^]^^^- 1 ]. (2.6) 

Now every [L] G PL(2, C) has a unique polar decomposition 

[L] = [U(A1 2 + \-t)\, (2.7) 

where ([U], A) G PP(2) x M 3 , A = V / TTA 2 , A = |A|, and Ti, r 2 , r 3 are the Pauli spin matrices. 
Hence Rati — PU{2) x IR 3 . Having chosen a basis {|r a } for su(2), we have a canonical 
identification PU(2) = SO(3) under which [U] is identified with the orthogonal transformation 
Adu : su(2) -> su(2). Hence Rati = 5*0(3) x 1R 3 , and the G action is 

(Pi, P 2 ) : (P, A) 1 — y (PiPP^ 1 , P 2 A). (2.8) 

From this we see that the G-orbits are level sets of A, generically diffeomorphic to SO(3) x S 2 
(when A > 0), the only exception being A = 0, which is diffeomorphic to 50(3). Physically, 
the lump corresponding to (P, A) G 50(3) x R 3 has maximum energy density at —A/A G 5 2 , 
sharpness proportional to A and internal orientation P. The A = lumps have uniform energy 
density. 

The following explicit characterization of G invariant kahler metrics on Rat x was established 
in [37] (see [3] for an alternative viewpoint, exploiting more directly the covering 5L(2, C) — > 
Rati): 

Proposition 1 Let 7 be a 50(3) x 50(3) invariant kahler metric on Rat x . Then 

7 = A x dA • dA + A 2 (\ ■ dA) 2 + A 3 a ■ a + A 4 (A • a) 2 + A t X ■ (<r x dA), (2.9) 

where A\,...,A4 are smooth functions of A only, all determined from the single function 
A\ = A(X) by the relations 

^2 = 7^ + ^, A 3 = \(1 + 2X 2 )A(X), A 4 = -1(1 + A 2 K(A). (2.10) 

Here o~i, o- 2 , o~ 3 are the left invariant one forms on 50(3) dual to the basis {|r a : a = 1, 2, 3} for 
su(2) = so (3), x and ■ denote the vector and scalar product on M 3 respectively and juxtaposition 
of one-forms denotes symmetrized tensor product. 

So symmetries and the kahler property determine the metric up to a single function A(X), 
which we may think of as the squared length of the vector d/dXi at the point (I3, (0, 0, A)) G 
50(3) x M 3 , corresponding to the rational map 

W{z)=vlz, A* = ^- (2.11) 
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For the L 2 metric, one finds that 



_ 327T/4// 4 - V log n-1] 

~ W^f • { } 

It follows from these formulae that (Rati, 7) has finite diameter and volume (16 3 7r 6 /6 accord- 
ing to Baptista [3]), is Ricci positive and has unbounded scalar and holomorphic sectional 
curvatures. Examining the large A behaviour of 7, one finds that the boundary at infinity of 
Rati is S 2 x S 2 . Geodesic flow in (Rat l5 7) was studied in detail in [36, 14], and turns out to be 
surprisingly complicated given the homogeneity and isotropy of physical space. In particular, 
lumps generically do not travel along great circles in S 2 . 

It is important to realize that any kahler metric on Rati invariant under the SO (3) x 50(3) 
action must have the structure of Proposition 1 for some function A(X). There is one other 
natural metric on Rati which we will have reason to consider. By identifying a rational map 
with the projective equivalence class of its coefficients, we may think of Rati as an open subset 
of CP 3 , 

Rati ^ CP 3 , aiZ + ° 2 ^ [ai, a 2 , a 3 , aA. (2.13) 

CZ3Z + 04 

We equip CP 3 with the Fubini-Study metric of constant holomorphic sectional curvature 4, 
whence Rati inherits a kahler metric, which we shall also call the Fubini-Study metric 7^5. 
It is not hard to show [37] that 7^5 is 50(3) x 50(3) invariant, and so has the structure of 
Proposition 1. One finds that the coefficient function is 

Afs = t+a? = Jn- i2M) 

Finally, given any 50(3) x 50(3) invariant kahler metric on Rati, it is convenient to define 
a second metric coefficient function, the squared length at ([I 2 ], (0, 0, A)) of 8 3 , where 9\, 9 2 , 3 
are the left-invariant vector fields on 50(3) dual to a%, 02, 03. This turns out to be 

B = ^(X 2 + A 2 )A + ^XA 2 A'. (2.15) 

For the L 2 and Fubini-Study metrics one finds, respectively, 

B = ^^K^ + ^bg/i-^ + l], 
11 

BfS ~ A(X 2 + A 2 ) 2 ~ A(2X 2 + 1) 2 - (2 ' 16) 



3 The Laplacian on Rat 
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We begin the computation of the Laplacian on functions on Rati by proving a lemma which 
generalizes the well-known expression for A in local coordinates. 
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Lemma 2 Let (M m ,g) be a Riemannian manifold of dimension m, {X{\ be a local frame on 
M, {z/j} be the associated coframe, and {/ij} be the associated basis for k m ~ l M , that is, 

Hx = v 2 Av 3 A---Av m , fi 2 = Vi A u 3 A ■ ■ ■ A v m , fi m = v x A v 2 A ■ ■ ■ A v m _ x . 

If all the (m — 1) -forms fa are closed, then the Laplacian on functions is 



\g\ 



i.j 



where g is the inverse metric and \g\ = det(g..). 

Proof: We use the summation convention on repeated indices, and define g % i = g{yi, Vj). The 
Laplacian on functions is A = — * d * d, where * denotes the Hodge isomorphism, so 



- * d(X t [/] *!/,) = -* d(XM} v[ 9 [ y g ij e ji2i3 ... im ^ A ■ ■ ■ A v im ) 
* X k [X^]^\g^]^^v k A v i2 A ■ ■ ■ A v im 



= - *X j [X i [f]y/\g\g t: >]u 1 Au 2 A---Au m , 

where we have used d(u i2 A ■ ■ ■ A u im ) = 0. Now the volume form is vol = \f\g\v\ A ■ ■ ■ A v n 
and *vol = 1, so 

Af=-X j [X i \f\y/\f\gW — 



\9\ 

as claimed. □ 

Note that any local coordinate basis X^ = d/dx l satisfies the conditions automatically, 
and that the formula for A reduces to the usual expression in this case. For our purposes it 
is convenient to use the (global) frame {d a ,9 a : a = 1,2,3} on Rati, where d a = d/d\ a and 
6 a are the left-invariant vector fields on SO (3) dual to o~ a . We thus require an expression for 
the inverse metric 7 relative to this frame. 

Proposition 3 Let 7 be a SO (3) x SO (3) invariant kahler metric on Rati determined, as in 
Proposition 1, by the function A(X). Then the inverse metric 7 is 

7 = dd-d + c 2 (\ ■ df + c 3 e-e + c 4 (a ■ of + c 5 x -(Oxd) (3.1) 

where C±, . . . , C5 are smooth functions of X alone satisfying 

X 2 + A 2 A 2 4 1 

Cl= A 2 A ' Cl + X2 ° 2 = 45 > C * = 1^A' C , 3 + A 2 C 4 = -, C 5 = C 3 . 

Here d = (8/6X1,8/8X2,0/8X3), = (0i,8 2 ,0 3 ), juxtaposition of vector fields denotes sym- 
metrized tensor product, and B is determined by A as in equation (2.15). 
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Proof: The same symmetry argument used to prove Proposition 1 (Proposition 3.1 of [37]) 
shows that every invariant symmetric (2, 0) tensor is of the above form, with C\, . . . , O5 some 
functions of A only. The formulae for C±, . . . , O5 result from explicit computation of ||dA a || 2 , 
||cx a || 2 and (d\ a ,o~b) a,b = 1,2,3, at the point ([I], (0, 0, A)) for general A, using the unitary 
frame {e a , Je a : a = 1, 2, 3} introduced in section 4.1 of [37]. For example, at ([I], (0, 0, A)), 

3 1 
IHI 2 = C 3 + A 2 4 = ^(a 3 (e a ) 2 + 3 (Je a ) 2 ) = 3 (Je 3 ) 2 = -. 

a=l 

□ 



Proposition 4 Let 7 be a SO (3) x 50(3) invariant kahler metric on Rati determined as in 
Proposition 1 by the function A(X). Then the Laplacian on (Rati, 7) is 



a/ = -^ I {e.ef + x.(dxe)f-^ 2 



A 2 A 

1 - 



AB 

1 8 f\ 2 A 3 A 2 df\ A 2 + A 2 



(A-0) 2 / 



A x d) ■ (A x d)f. (3.2) 



AX 2 A 2 B d\\ 4 dXJ X 2 A 2 A 

Proof: We work with the frame {d a , 8 a : a = 1, 2, 3}. Note that d(dA a ) = and 

dcr x = cr 2 A 03, do- 2 = 03 A 0i, d0 3 = 01 A 2 , 

so any wedge product of five of these one forms is closed. Hence this frame satisfies the 
conditions of Lemma 2 which, with Proposition 3, immediately gives 

a/ = -^^d-(y^\(c 1 df + c 2 x(\-d[f}))) + d-(- v ^\^\xe[f}) 

+0 ■ (vl^yA x d[f]) + 6 ■ (VW\(C 3 0[f] + 4 A(A • ©[/])))}. (3.3) 

Now the volume form on (Rati, 7) is [37] 

vol = -A 2 BdX 1 A dA 2 A dA 3 A o x A a 2 A 03, (3.4) 

whence one sees that \/[y\ = ^A 2 B for this frame. Substituting this, and the formulae for 
Oi, ...,05 (Proposition 3) into equation (3.3) yields, after some straightforward manipulation, 
the formula claimed. □ 

As a check on our formula, we should verify that it is consistent with the 50(3) x 50(3) 
symmetry of 7. That is, the operator A must commute with all Killing vector fields on 
(Rati, 7). There are six independent Killing vector fields on (Rat 1; 7), three generating the 
left 50(3) action and three generating the right 50(3) action. Recall that the left action 
on Rati - 50(3) x 1R 3 acts by left translation on 50(3) and acts trivially on M 3 . Now left 



S 



translations on a Lie group are generated by right invariant vector fields. So, let £ a , a = 1, 2, 3, 
be the right invariant vector fields on 50(3) coinciding at I with 9 a . Since the left and right 
actions commute, [0 a ,6>] = f° r & U a i b so we see immediately that A as in formula (3.2) 
commutes with each £„. The right action of 50(3) on Rati — 50(3) x M 3 acts by right 
translation on 50(3) and by the fundamental representation on IR 3 . Hence the Killing vector 
fields generating this action are X a = 9 a + <3> a where $ a = e a b c Xbd c . Note that {O a }, and 
{X a } each satisfy the angular momentum algebra, 

[0a,6b] = -eabcdc, [&a, &b) = -^abc^c, [X a , X b ] = -€ abc X c , (3.5) 

and that [9 a , $&] = 0, so 

[X a , X-X] = [X a , 0-0] = [X a , $ • *] = 0. (3.6) 

Clearly all these vector fields annihilate functions of A and commute with d/d\. It follows 
that X a commutes with the first, fourth and fifth terms of (3.2). To deal with the second 
term, note that 

[X a ,X- (d x 0)] = [X a ,0- *] = l -[X a) X- X-00-$$] = O. (3.7) 
Finally, X a commutes with the third term since 

[X a , A • 0] = [9 a + CabcXbdc, Xd9d] = Xd{ — € a de9 e ) + ^abd^b9d = 0. (3.8) 

4 Casimir energy 

The dominant term in Moss and Shiiki's Hamiltonian arising from the degrees of freedom 
orthogonal to the moduli space is the Casimir energy c io , which is formally given by 

i 

where are the frequencies of oscillation of the normal modes of the static solution. In the 
field theory context, the above sum is infinite and divergent, so we must regularize it in some 
way. We return to this issue below. First we set about computing the frequencies cuj in the 
case of interest, a single CP 1 lump on 5 2 . To achieve this, we will make use of some standard 
results in the stability theory of harmonic maps, which we begin by briefly reviewing. This 
material is treated in detail in [40, ch. 5]. 

4.1 The spectrum of the Jacobi operator 

Recall that given a map tp : (M, g) — » (N, h) between Riemannian manifolds, its Dirichlet 
energy is 

E = \j (4.2) 
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and the map is harmonic if it is a critical point of E. This is directly relevant to us since, with 
the choice (M,g) = (N, h) = the round sphere of radius 1, E coincides precisely with V, the 
CP 1 model's potential energy functional. Hence, CP 1 lumps are harmonic maps S 2 — > S 2 . 
Given a harmonic map tp : M — > X, one defines its Hessian, a symmetric bilinear form on 
r(<^ -1 TX) (the space of smooth sections of the vector bundle over M whose fibre over p £ M 
is the tangent space T V ^N), as follows. Let (p Stt '■ M — > X be a smooth two-parameter 
variation of (p (so ^ ,o = V?), with d s ip s j\ s=t=0 = X, d t <f s ,t\s=t=o = Y £ F^^TN). Then 

(4.3) 

s=t=0 

One says that ip is stable if Hess^X, X) > for all X. Clearly, if tp minimizes E in its 
homotopy class, as in our case, it is stable. 

Associated with Hess^, there is a self-adjoint, elliptic, linear differential operator J;^ : 
r(<£> -1 TiV) — > r(y? _1 TX), which is known as the Jacobi operator and is defined by 

Hess v (X, Y) = (X, jf v Y) L * = [ h(X, f^Y). (4.4) 

J M 

If M is compact, the harmonic map <p is stable if and only if the spectrum of J?^ is non- 
negative, and this spectrum is discrete, each eigenvalue having finite multiplicity. Any map 
%p : M — > N which is sufficiently close pointwise to a harmonic map (p can be uniquely written 
ip = exp^X, where ||X|| is pointwise small, and then 

Ety) = E{<p) + ^Hess v (X, X) + 0(X 3 ) = E{<p) + ±(X, f^X) L , + 0(X 3 ), (4.5) 

whence it is clear that the eigenvalues of are precisely u 2 , the squared frequencies we 
require to compute . 

There is an explicit formula for the Jacobi operator of a general harmonic map (p : (M, g) — > 

(N,h), 

= A v - & v , (4.6) 

where is the rough Laplacian on T(ip~ 1 TN), and is a certain section of End(ip~ 1 TN) 
constructed from the curvature tensor R N on X. Explicitly, given a choice of local orthonormal 
frame Ei, E 2 , . . . , E m on M, 

A,Y = -trV-V-Y = -f^hl^lY) -VI Me y\ , 

8=1 ^ 1 ) 

III 
8=1 

where V M , V N are the Levi-Civita connexions of M,N respectively, and is the pullback 
to ip>~ x TN of VX In the case of interest to us, is somewhat easier to handle than A^, 
owing to the following proposition: 
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Hess^(X,y) 



d 2 



ds dt 



(4.7) 
(4.8) 



Proposition 5 Let ip : (M n ,g) — > (N n ,h) be a weakly conformal mapping between Rieman- 
nian manifolds of equal dimension, and (N n , h) be Einstein with scalar curvature k (necessarily 
constant). Then 

Hsu, 



n z 

where £ = ^\\dip\\ 2 G C°°(M) is the Dirichlet energy density of ip. 

Proof: Recall that <p is weakly conformal if there exists a smooth function / : M — > R such 
that h(dtpX, d<pY) = f 2 g{X, Y) for all vector fields X, Y on M. Let E x , . . . , E n be a local 
orthonormal frame on M. At all points p G M where f(p) = 0, dip = so the desired equality 
holds trivially (both and £ vanish). At all other points, || d^.E'^ || 2 = -£ > independent 
of % and h(dtpEi, d<pEj) = for i ^ j, so 

/ n 

Ei = J—dipEi, i = l,2,...,n 
form an orthonormal basis for T^^N. Hence, for all X, Y G T^^N, 

h(X, ^Y) = £ h (X, R N (Y, d<pEi)d(pEi) = —J2 M*, rN ( Y > E i) E i) 

i=l i=l 

= -—f2h(E l ,R N (Y,E l )X) = —R\c(Y,X) = —-h(X,Y) (4.9) 
n z — ' n n n 

i=l 

by a standard symmetry of R N [41, p. 58] and the Einstein property of (N, h). □ 

Note that every ip G Rati is holomorphic, hence conformal, and the unit two-sphere is Einstein 
with k — 2. Hence 

& v = £Id (4.10) 

in this case. Thinking of Jfi^ as a quantum Hamiltonian acting on "wavef unctions" on the 
two-sphere (sections of (p~ x TN, really), the effect of the curvature term is to add a potential 
well equal to minus the classical energy density of the lump. What is this energy density? By 
50(3) x 50(3) invariance, the spectrum of J?^ can depend only on A, so it suffices to consider 
only the one-parameter family of rational maps 

ip : z ^rW = fi(X)z = ~^—^ z ( 4 - n ) 

corresponding to the curve ([I2], (0, 0, A)), A > 0, in Pf/(2)xR 3 . It is convenient to parametrize 
these by n G [1, 00) rather than A G [0, 00). In terms of the usual polar coordinates on both 
domain and codomain spheres, the map (4.11) is 



¥>:(M)H>(/ M (0),0), 



where / M (0) = 2 cot" 1 ^cot Q . (4.12) 



The pair E\ = d/d8, E 2 = cosec 6d/d<f) is orthonormal on S 2 , whence, by conformality, one 
sees that 

£=\\d V E 2 \\ 2 = S ^J^l. (4.13) 
sin (7 
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We turn now to the explicit computation of the rough Laplacian A^,, for ip of the form 
(4.11). Note that E\ = Eiotp, E 2 = E 2 otp gives an orthonormal pair of sections of p~ l TN . Any 
section Y of (p~ 1 TN can be uniquely written Y = Yi(9, <f))Ei + Y2(9, (j))E 2 . We seek an expres- 
sion for A v as a differential operator acting on the pair of smooth functions (Yi, Y 2 ). In fact, 
we can deduce all we need once we know how acts on sections of the form a(9) cosm<j)Ei 
where m G N. A straightforward but lengthy computation, presented in the appendix, shows 
that 

A (fi (a(9) cos m(pEi) = (@ m a) cos m$E\ + Q m a sin m$E 2 (4-14) 

J2 J ^,2 i _2 



where -lm ~ cot 9 Ta + m + T Q U 

d9 2 d9 sin 9 



Q m = 2m^% (4.16) 
sin 9 

Since N is kahler, V w commutes with J N , the almost complex structure on N, namely 
[V^, J N ] = 0, and hence [A^, J N ] = also. So we conclude immediately that 

A ^(a(9) cos m(pE 2 ) = A„( J N a(9) cosm^Ej = J N A ip (a(9) cosm^) 

= — Q m asmm(f)Ei + (@ m a) cosm(j>E 2 . (4-17) 

Consider now the 5*0(2) action on C°°(S 2 , S 2 ) given by 




p (->• (p a = R a o ip o R_ a where R a = I sin a cos a 1, ael. (4.18 



Note that E(ip a ) = E(ip) and each of the maps (4.11) is fixed under this SO (2) action. It 
follows that and hence, by Proposition 5, A^, commute with the induced action on 
Tip-'TN), 

Y H> R a (Y) = dR a o Y o R_ w (4.19) 
or, in terms of polar coordinates, 

R a (Y 1 (9, <f>)E x + Y 2 (9, <P)E 2 ) = Y 1 (9, - a)E, + Y 2 {9, - a)E 2 . (4.20) 

Hence 

A ip (a(9) smm(pE 2 ) = A ip (R^a(9) cosm<j)E 2 ) = R_zi_ A ip {a{9) cosm<j)E 2 ) 
= Rx (—Q m a sin m(j)Ei + (@ m a) cos m(f)E 2 ) 

= Q m a cosm(f)Ei + (@ m a) sinm(pE 2 , (4-21) 

and similarly 

A v (a(9) smm(j)Ei) = (@ m a) sinmipEi — Q m acosm<j)E 2 . (4.22) 
To summarize, J?^ leaves each infinite dimensional subspace I m C T(p~ l TN), m G N, 

Im = {a±(9) cos m(j)Ei+a 2 (9) sin m(j)E 2 +as(9) cos m(j)E 2 + ai(9) sin m(pEi : a G C°°((0, vr), M 4 )} 

(4.23) 
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invariant, these spaces span T(if 1 TN), and on /„ 



a 2 
a 3 
\ a 4 / 



/ 



V 



Qm 






£ Qn 

-&m 






£ 






~Qr 



£ 






-Qr 



£ J 



a 2 
\ a 4 / 



(4.24) 



This diagonalizes after a simple change of coordinates. Let a\ 
«3 = a>i — d2 and a 4 = 0,3 — a 4 . Then 



ai + a 2 , «2 = «3 + «4, 



( «i \ 

\ «4 / 



/ 



V 



Note that 



Qr 



o 





and Q-ri 



\ / «i \ 

+ Q m -£ a 2 

o ^ m -g m -£ a 3 

3> m - Q m - £ J \ a 4 ) 

(4.25) 

—Qm, so spec is the union of the spectra of the 



Sturm-Liouville operators 

S m Qr 

where 

U m (9) - 



£ 



d 2 



d6 2 



cote 



d 

1e 



U m (9), 



m G Z, 



m 



l + 2cos 2 f^9) + 2mcosf fl (9) 
sin 2 9 



(4.26) 



(4.27) 



and each eigenvalue occurs with double multiplicity. In terms of physics, we may think of this 
as the energy spectrum for "spin 0" (i.e. independent) states of a point particle moving on 
S 2 in the 5*0(2) invariant potential well U m {9). The spectrum of each S m may be computed 
numerically using the shooting method described in [14]. 

There is one value of fi for which spec may be computed exactly, namely \i = 1. 
Here things simplify considerably, because the corresponding rational map is the identity 
map S 2 — > S 2 . The Jacobi operator for the identity map on a general Riemannian manifold 
(M n , g) was studied in detail by Smith [35]. The key simplification is that one has a canonical 
identification Id^TiV = T*M, obtained by identifying the section Y of Id _1 TiV = TM with 
the one form b Y = g(Y, •). This is useful because there is a Weitzenbock formula relating the 
rough Laplacian Ay to the Hodge Laplacian A on one-forms [40, p. 161] 



Ai d Y = (A b F - Ric(F, •)). 
In the case where M n is Einstein, with (constant) scalar curvature k, 



Ric(F,-) 



Y. 



n 



(4.28) 



(4.29) 



which together with Proposition 5 and the observation that £ — §, constant, for Id, gives the 
following formula for J?\& on an Einstein manifold, originally due to Smith, 



id 



A , 2K 

Ab . 

n 



(4.30) 
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In our case, n = k = 2, so = fj A b — 2, that is, the spectrum of J?^ in the special case 
/i = 1 is the spectrum of the Hodge Laplacian on one-forms on S 2 , shifted down by 2. This, 
in turn, can be related to the spectrum of the Laplacian on functions on S 2 , as follows. 

Let Q p denote the space of smooth p- forms on S 2 , and 5 : Q p — > Q p ~± denote the coderiva- 
tive, so that the Hodge Laplacian is A = 5d + d5. Every one-form on S 2 has a unique Hodge 
decomposition 

Yi = dYo + SY 2 (4.31) 

into exact and coexact components (the harmonic component vanishes as H l (S 2 ) = 0). Now 
df2 and 5Q 2 are L 2 orthogonal subspaces of Q\ and [A,d] = [A, 5] = 0, so the spectral 
problem for An x decomposes into two sub-problems, for A|df2 and A\sn 2 - Since Y\ is coexact 
if and only if *Yi is exact, and * : Qi — > fii is a L 2 isometry, the spectra of A|dn an d ^\sn 2 
coincide. Hence spec is spec A|dn with double multiplicity. In fact, 

spec A| dno = spec A no \{0}. (4.32) 

To see this, let v £ spec Aq and Y be the corresponding eigenfunction. Then 

A dYo = dAYo = dz/Y = ^dY (4.33) 

so if dY 7^ (that is, v ^ 0), then v e spec A|dn . Conversely, let v e spec A|dc an d dYo be 
the corresponding eigenform. Then v ^ (since there is no harmonic one- form on S 2 ), so 

= AdY - udY = d(AY - uY ) 
AY — uYq = c, constant 

AYJ-^YJ = 0, where Y ' = Y + £ (4.34) 

and hence v £ spec Aq . The spectrum of Aq is well known. 
We conclude that, at /x = 1 (equivalently A = 0), 

spec = {£(£ + 1) - 2 : £ £ Z + } 
multiplicity {£{£ + 1) - 2) = 4£ + 2. (4.35) 

Since we know the eigenvalues (and eigensections) of at A = 0, we can use these as seed data 
for the numerical shooting method, starting at A = and increasing A in small steps. In this 
way we can numerically construct curves oj 2 {\) showing how the different eigenvalues vary with 
lump sharpness A. Figure 1 shows such curves for the lowest 48 eigenvalues. Recall that the 
eigenvalue has multiplicity 6 (the dimension of the moduli space), and that other eigenvalues 
always have multiplicity (at least) 2 due to the symmetry under J N (or, equivalently, due to 
the block structure of J?^ on I m ). 

It is interesting to examine the A — > oo behaviour of u) 2 (\). The pointwise limit of <p : 
S 2 ->• S 2 as A ->■ oo is 

«> (p) = { (0 ' ' 1} ^t ' '" 1 ) (4 36) 

^°° [P) \ (0,0,-1) p = (0,0,-1) 1 j 

that is, if oo is constant almost everywhere. A sensible guess for the limiting spectrum would, 
therefore, be the spectrum of the Jacobi operator at a constant map, which is known to coincide 



14 











5 



10 



15 



X 



Figure 1: The dependence of the eigenvalues of the Jacobi operator Jp^ on the lump sharpness 
A. Note that the eigenvalues interpolate between £{£ + 1) — 2, £ G Z + , at A = and £(£ + 1), 
£ e N as A ->> oo. 

with n copies of the spectrum of the Laplacian on functions, where n is the dimension of the 
codomain [40, p!60]. In this case 



The numerics suggest that this guess for the limiting spectrum is very nearly correct. Specif- 
ically, the eigenvalues do tend, as A — > oo to eigenvalues of ^const, and, apart from the 
eigenvalues and 2, those eigenvalues tending to £(£ + 1) have total multiplicity 4£ + 2. So 
the guess is wrong only in that it predicts the multiplicity of the limiting eigenvalue to be 
2 rather than 6 (as it must be given the dimension of Rati) and the multiplicity of 2 to be 6 
rather than 4 (as found numerically). 

4.2 The regularized Casimir energy 

Having computed spec Jf^, we must now try to make sense of the Casimir energy ^(A) in (4.1) 
which, as it stands, is divergent. It is conventional to reset the zero of potential energy so that 
the total zero-point energy of the vacuum is 0. For our purposes, y?oo is (almost everywhere) 
the vacuum, so it is convenient to define 



Spec j£ 'const 

multiplicity (£(£ + 1)) 



{£{£ + !) : £eN} 
4£ + 2. 



(4.37) 




(4.38) 



8=1 
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where the eigenvalues 0Ji(X) 2 > are arranged in such a way that Ui(0) 2 is nondecreasing. 

One option would be to renormalize the spectrum (4.38) numerically, using for example 
the heat kernel approach. Here, we will employ a different semi-analytical approach. We first 
discuss the finite sums ^ for special values of k and then describe how we regularize the 
diverging sum as k — > oo. We choose k — 10, k — 24 and k = 42, which include the lowest 1, 2 
and 3 eigenvalues of respectively. Note that, since we have eigenvalue crossings in figure 1, 
this amounts to the lowest k normal modes at A = 0, but not at large A (where the eigenvalue 
ordering has changed). If we were to define ^(A) as the sum of the frequencies of the lowest 
k normal modes at each A, the function ^ would not be smooth (it would have corners where 
eigenvalues cross). In effect, we are making a large but finite-dimensional approximation to 
the quantum field theory, in which the wavefunction is a function on a vector bundle over 
Rati, whose fibre over ip is a union of low- lying eigenspaces of a ^' lp . We are choosing these 
eigenspaces so that they vary smoothly over Rati. The price for this is that they are, towards 
the boundary of Rati, not quite the lowest energy eigenspaces available up to dimension k. 




5 10 15 

X 



Figure 2: The one-lump Casimir energy as a function of lump sharpness A for three different 
cut-offs: solid curve ^io(A), dashed curve 1^24 (A), dotted curve ^^42 (A), where ^ is defined in 
(4.38). The values k — 10, k — 24 and k = 42 correspond to the lowest 1, 2, and 3 eivenvalues 
of Jia-, respectively. 

Plots of ^io, ^24 and ^ 42 are presented in figure 2. The curves for ^ 2 4 and ^42 have been 
rescaled vertically, by a factor of \ and | respectively, to make comparison with ^10 easier. 
Note that these three functions are qualitatively very similar. In fact, the Casimir curves ^10, 
^24 and ^42 are approximately self-similar up to a factor which diverges as k — > 00. We have 
opted to use the self-similarity and regularize the diverging factor (the depth of the well at 
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A = 0). So we define the approximate renormalized Casimir energy ^(A) to be 



where C* is the renormalized Casimir energy of the A = lump. This can be computed exactly 
using zeta function regularization, because, as we have seen the spectra of the Jacobi operator 
for the identify map (A = lump) and the constant map (the vacuum) are known exactly. 



4.3 Zeta function regularization 

For A = 0, the spectrum of the Jacobi operator is known explicitly, namely Uq = with 
multiplicity fio = 6 and uf = 1(1 + 3) with multiplicity \ii = 4/ + 6 for / = 1, 2, . . . This enables 
us to calculate the Casimir energy using zeta function regularization. The key idea is to write 
down the corresponding zeta function 

oo 

C( S ) = $>Kr , (4.40) 
1=1 

leaving out the zero modes. The zeta function (4.40) is absolutely convergent as long as the 
real part of s is sufficiently large. In this region, ((s) can be viewed as an analytic function 
of s. We are interested in the Casimir energy which corresponds to the value s — — |. In this 
case, the formal sum in (4.40) is divergent, however, ((s) defined as analytic continuation is 
well-defined. It is convenient to rewrite ((s) as the following sum: 

oo oo 

cc«) = £ 21(1(1 + 3)r + E 2 c + 3 W + 3 ))" s ( 4 - 41 ) 
i=i i=i 

then we can use formula (5.8) in [9, p. 122] to obtain 

C(-|) = -2.373. 

The corresponding calculation for the vacuum, whose spectrum is uf = 1(1+1) with multiplicity 
Hi = 41 + 2, leads to the zeta function 

oo oo 

COO = $>(/(/ + 1))-* + £)2(J + l)(l(l + I))"* (4-42) 
1=1 1=1 

and yields 

CvacH) = -0.530. 

See also equation (5.34) in [9, p. 126]. Hence, the Casimir energy of the A = lump on the 
unit two-sphere can be evaluated using 

\ (c(-|) - CvacH)) = -°- 921 - 
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Hence 

C* = |^(0)| = 0.921. 

This is the energy scale we used to set the renormalized energy scale for our numerically 
generated Casimir energy function ^(A) in (4.39). 

Our approximation to the Casimir energy ^(A) is now given by the rescaled and shifted 
curve ^io (A) in figure 2. It is non-singular and appears to be smooth. It is worth summarizing 
the approximations involved. Our calculation relies on a conjectured self-similarity of the 
curves in figure 2. We also assume that we can neglect the effects of crossing modes. In 
particular, we assume that the zeta-function regularization can be performed pointwise when 
we regularize the Casimir energy at A = and A = oo. Furthermore, we assume that the 
^(A) — >• as A because the lumps converge almost everywhere to the vacuum in that 
limit. 

The obvious alternative to our approach, namely a numerical evaluation using heat kernel 
or zeta function regularization as in [29, 27], is beyond the scope of this paper. 



5 The energy spectra 



In the following, we describe how to calculate the spectrum of the Laplacian (3.2). In order 
to make use of the physics literature on this topic we rewrite the Laplacian using angular 
momentum operators. Note that the operators J = — i6, L = —iX x d, T = J + L and 
K = —i£ all satisfy the canonical commutation relations for angular momenta, namely, 



iX 



3- 



X = J, K, L, T 



and cyclic permutations. Recall that generates the left 5*0(3) action, on the target 5 2 , 
and iT generates the right 50(3) action, on the physical 5 2 , so we refer to these operators 
as isospin and spin respectively. Making use of these operators, one sees that the Laplacian 
(3.2) is 



d fX 2 A 3 A 2 diP 



AX 2 A 2 B dX 
2 



4 



A 2 A 



J z iP + TV 



dX 
2 

A^ 



+ 



X 2 A 2 A 
A 2 A 



L 2 i/j 



45 



A • T) V 



Here, we replaced the L • J term using the convenient operator identity 

2L J = T 2 — L 2 — J 2 . 



(5.1; 



(5.2) 



We have already shown that A commutes with T and K, as it must, by 5*0(3) x 50(3) 
invariance. In fact, the operators A, |T| 2 , T 3 , |K| 2 , K 3 , P, where P : A i— > —X is the parity 
operator, all mutually commute, so we can seek simultaneous eigenstates of these six operators. 
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5.1 Spin-isospin interchange symmetry 

Let t(t+l) and k(k + 1) be the eigenvalues of |T| 2 and |K| 2 . It is a somewhat surprising fact 
that the spectrum of A is invariant under interchange of t and k. Even more surprising is 
that this remains true even if we give the physical and target spheres different radii. The key 
to seeing this is the following isometry of any 50(3) x 50(3) invariant kahler metric on Rat^ 

Proposition 6 Identify Rati = 50(3) x R 3 with T50(3) via (R, A) = (\-0)(R) E T R S0(3). 
The mapping f : 50(3) — > 50(3), f(R) = R^ 1 induces a mapping df : T50(3) — > T50(3). 
Then df is an isometry of any 50(3) x 50(3) invariant kahler metric on TS0(3). 

Proof: The identification Rati = T50(3) amounts to thinking of A as a vector in the Lie 
algebra so (3) (skew 3x3 real matrices) and identifying Tr50(3) with so (3) = T e 50(3) by 
left translation. Hence, the map in question is 

df:(R } X)^(R-\-Ad R \). 

Let q(t) = (R(t),X(t)) be a curve in Rati with q(0) = (R(0)fl,v), £l,v G so(3). Then, with 
respect to the 50(3) x 50(3) invariant kahler metric determined by the functions Aj(A), as 
in Proposition 1, 

||g(0)|| 2 = A 1 \v\ 2 + A 2 (X ■ vf + Az\n\ 2 + A 4 {\ ■ SI) 2 + A x v ■ (A x ft). 

The image of this curve under df is q(t) = (i? _1 (t), — i?(t)A(t)i? _1 (t)), which has g(0) = 
(R-\0)n,v), where 

(n,v) = -(Ad R n,Ad R (v - [A, n]„( 3 )))- 

Recalling that the Lie bracket on so (3) coincides with minus the vector product on M 3 under 
the natural identification so(3) = M 3 , and the adjoint action of 50(3) on so(3) coincides with 
the fundamental action on IR 3 , we find that 

||g(0)|| 2 = + A x Q\ 2 + A 2 {\ ■ v) 2 + A 3 \Q\ 2 + A 4 (A • fl) 2 - A x {v + A x Q) ■ (A x n). 

= H(o)f. 

Hence df is an isometry, as claimed. □ 

Proposition 7 Given Eel and k, t e N, denote by XE,k,t the simultaneous eigenspace of 

H=±A + V, |K| 2 , |T| 2 

with eigenvalues E, k(k + 1) and t(t+ 1), where V : Rati — > K is any 50(3) x 50(3) invariant 
potential. Then XE,k,t an d XE,t,k have equal dimension. 

Proof: Under the identification Rati = TS0(3), the left and right actions of 50(3) on Rati 
coincide with the natural left and right actions of 50(3) on T50(3). Denote by df* the induced 
map on L 2 (T50(3), C), df *ijj = ip o df . Then, since df is an isometry and preserves the length 
of A, Hdf* = df*H. Furthermore, df interchanges the left and right 50(3) actions on T50(3), 
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so Kdf* = -df*T and Tdf* = -df*K. Hence, df* : X E ,k,t -> X E , t)S; and d/* : X E)t , fe -> 

for all E,k,t. Now rf/ 2 = Id, so (df*) 2 = Id, and hence df* is invertible, so the eigenspaces 

XE,k,t and X E ,t,k are isomorphic. □ 



Note that this k,t interchange symmetry holds for any invariant kahler metric on Rati, 
so is not special to the case where the domain and target two-spheres have equal radius. 
Also note that it relies heavily on the kahler property of the metric, and does not follow from 
5*0(3) x 50 (3) invariance alone (or, indeed, just invariance and hermiticity). This is a "hidden" 
symmetry which one cannot see directly from the field theory. It is an interesting question 
whether a direct physical argument can be given to explain this symmetry. A construction 
involving a supersymmetric extension of the model would be a natural candidate. 

5.2 Reduction to a Sturm-Liouville problem 

The wavefunction is a function ip : SO (3) x M 3 — > C. According to the Peter- Weyl theo- 
rem [16], the matrix elements of all irreducible unitary representations of the group 5*0(3) 
form an orthonormal basis for L 2 (50(3), C). Recall that such representations are labelled 
by k G {0,1,2,...}, and that the matrix elements of the k representation are functions 
7Tj^ kg : 50(3) — > C where —k < js,k 3 < k. We have chosen the symbol k to label the 
representations because 

l K l^l = UNl = ^ + l)-jl- (5-3) 
Further, we may choose the basis for C 2fc+1 so that 

IT ( k ) 7 ( k ) T ( fc ) ' ( fc ) / r A\ 

which is why we have labelled the matrix elements with j'3, k^. For the M 3 dependence of ip we 
use spherical polar coordinates and expand the angular dependence in spherical harmonics. 
Hence we may express the wavefunction as follows, 

00 00 k k I 

^ = EEE E E a% k3>l3 (\)\k,l, j3 ,k 3 ,l 3 ), (5.5) 

fc=0 1=0 j 3 =-k k 3 =-k h=-l 

where 

\k, l,j 3 , h, h) = x^Yi^, (5.6) 

and Yy 3 : 5 2 — > C denote spherical harmonics. 

Recall that H, |K| 2 ,i^3, |T| 2 ,T3 are mutually commuting, so we may solve the eigenvalue 
problem for H on each simultaneous eigenspace of |K| 2 ,i^3, |T| 2 ,T3 separately. Clearly, H 
is independent of k% and £3 so we may, and henceforth will, without loss of generality, set 
^3 = £3 = 0, remembering to multiply all degeneracies by (2k + 1)(2£ + 1) to account for the 
other values of (k 3 ,t 3 ). Furthermore, by Proposition 7 we may, without loss of generality, 
assume that k < t, doubling the multiplicity if k < t. Recall that T = J + L and that J and 
L satisfy the angular momentum algebra. Hence a basis {\k,t,l) : \k — t\ < / < k + t} for 
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the (k, 0, t, 0) eigenspace can be constructed, using Clebsch- Gordon coefficients, on which the 
operators |J| 2 = |K| 2 , |T| 2 , |L| 2 act naturally 



\3\ 2 \k,t,l) = k(k + l)\k,t,l) (5.7) 
|T| 2 |M,0 = t(t + l)\k,t,l) (5.8) 
|L| 2 |M,Z) = l(l + l)\k,t,l), (5.9) 



see Appendix B for details. 

By expanding the wavefunction as ip = Yl\=\k-t\ a iW \^,t,l), we may express the Hamil- 
tonian H = | A + F as a m x m matrix of differential operators acting on the vector function 
a : [0, oo) — > C m , where m = 2mm(k,t) + 1. Its structure is 

Ha=- Pl (\)^ (^ P2 (\)^+p 3 (\)M 1 a+p 4 (\)(k(k+l)+t(t+l))a+ P5 (\)M 2 a+V(X)a (5.10) 

where Pi(\) are rather complicated but explicitly known functions of A (see Proposition 4) 
and Mi, M 2 are m x m matrices corresponding to L 2 and (A • T) 2 respectively. By our choice 
of basis, Mi is diagonal and has entries /(/ + 1) with / running from \k — t\ (top left) to k + 1 
(bottom right). The matrix M 2 corresponding to operator (A • T) 2 is the only non-diagonal 
term in the Hamiltonian and is discussed in more detail in appendix B. It mixes states of 
different I, which differ by 2. This leads to a natural chess board structure. Hence, by a 
reordering of the basis vectors, the operator can be written in block diagonal form with one 
block corresponding to the states with even I and the other block to odd I. These blocks 
correspond to the decompostion of the (k, 0, t, 0) eigenspace into P = +1 (even /) and P = — 1 
(odd /) parity eigenspaces noting that PY^ m = (— l) l Yi im . In summary, for fixed (k, t), the 
eigenvalue problem for any Hamiltonian of the form H = |A + V(X) reduces to a matrix- 
valued Sturm-Liouville problem on [0, oo) of dimension m = 2min(/c,t) + 1. 



5.3 Boundary conditions 

We can now address the eigenvalue equations 

H ia = Ea, (5.11) 

where E is the energy eigenvalue, and the Hamiltonians Ho, Hi and H 2 are given by (1.3). As 
just explained, the spectral problem reduces to a sequence of matrix-valued Sturm-Liouville 
problems indexed by t G {0, 1, 2, . . .} and k G {0, 1, . . . , t}, where each subproblem has dimen- 
sion 2 min(/c, t) + 1. In order to calculate the spectrum of the Hamiltonians H , H 1 and H 2 we 
have to derive not only the relevant differential equations but also the appropriate boundary 
conditions. For many important examples in mathematical physics, the boundary conditions 
are determined solely by the requirement that the wavefunction be L 2 finite. However, in 
our case, L 2 finiteness is not always sufficient, and we have to apply the theory of singular 
Sturm-Liouville equations, following [30]. 

The asymptotic behaviour of the L 2 metric on Rati can be obtained by direct calculation: 

A ( x ) ~ IT asA-^0, A(X) ~ w asA^oo, 

£(A)~f asA^O, B(\)~^zA as A oo. 1 > 
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Using the above limits, we obtain the leading order equation for (Hq — E)f = as X — > oo, 



d 2 f , Idf 4 fi y \ 

T^T-T^ AT / = 0» wlth A =T> ( 5 - 13 ) 



<9A 2 A<9A A 2 V ) J ' A 
and the leading order equation for (Hq - £)/ = as A 

d 2 f , 2df 1(1 + 1) 



dX 2 XdX A 2 



/ = 0. (5.14) 



For A equation (5.11) for Hq has a regular singular point. The asymptotic form of the 
solutions can be derived from (5.14) and is given by 1 

f(X) = c 1 X l + c 2 X- 1 - 1 . (5.15) 

Note that the leading order asymptotic behaviour is independent of the energy eigenvalue E. 
We are interested in solutions that are L 2 finite which leads to the condition 

/ 2 vol < oo. (5.16) 

Using polar coordinates for (3.4) we obtain the two asymptotic behaviours, 

/ 2 vol ~ ~ Cl X 2l+2 and f\o\~~c 2 X- 2 \ 

for suitable constants c\ and c 2 . Hence I > corresponds to the limit-point case, since only the 
first solution is L 2 finite, see e.g. [30] for further details. However, for / = the situation is 
slightly more subtle. Both asymptotic solutions are L 2 finite. This is known as the limit-circle 
case. 2 The Hamiltonians H 1 ,H 2 lead to the same asymptotic behaviour since the curvature 
and Casimir energy are finite as A — > 0. 

For A — > 00 it is convenient to analyse the boundary conditions in a different basis so that 
the operator (A • T) is diagonal, 

(a-t) 2 / = p 2 /, 

where the integer k satisfies — min(A;, t) < p < min(A;, t), see appendix B. Then the asymptotic 
solution follows from (5.13) and is given by 

/(A)= Cl A" 2p + c 2 A 2p , (5.17) 

for p 7^ and 

/(A) = Cl + c 2 log(A), (5.18) 

for p = 0. Again the leading order behaviour is independent of the energy eigenvalue E. 

Both Hi and H 2 contain the scalar curvature function k(X) which is known to diverge to 
infinity as A — > 00. In fact, using the formula obtained in [37], and accounting for the change 



1 Higher order terms in the expansions for A — > and A — > 00 have been calculated in [24] for Hq with k = 0. 
2 The same subtlety occurs when solving the Schrodinger equation for the hydrogen atom in spherical polar 
coordinates. 
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in normalizations (recall we are giving the domain and target spheres unit radius), we find the 
asymptotic formula 

1 A 4 

re (A) ~ as A -»• oo. (5.19) 

v ; 167rlogA v ; 

This adds a term to (5.13), namely, 



dX 2 XdX X 2 V V / 21og(A) 
For p ^ we can solve (5.20) in terms of modified Bessel functions 



/(A) = dVlogA^^ (2plog(A)) + c 2V 4oiA/^ /2 (2plog(A)) . (5.21) 

With the help of asymptotic expansions, see e.g. [1], it can be shown that the leading order 
is identical to (5.17). For p = we obtain 

/(A) = c 1 log(A)^5v / 3 + C2 i og (A)l+^A (5.22) 

Since the Casimir energy is bounded, the asymptotic behaviour of H 2 as A — > oo is the same 
as Hi. 

As A — ?• oo the L 2 finite condition (5.16) again leads to two asymptotic behaviours, 

/ 2 vol ~ ci log(A)A~ 2p - 5 and / 2 vol ~ c 2 log(A)A 2p - 5 , 

for suitable constants c\ and c 2 in the case p > 0. For p — 0, the curvature term has an 
influence on the asymptotic behaviour. For Hq, we obtain 

/ 2 vol ~ ci log(A)A~ 5 and / 2 vol ~ c 2 log(A) 2 A~ 5 , 



whereas for Hi and H 2 the asymptotic behaviour is 

/ 2 vol ~ ci log(A)^~^A" 5 and / 2 vol ~ c 2 log(A)5 + ^A -5 . 

Hence for p < 2 both solutions lead to L 2 finite solutions, resulting in a limit-circle case. 

In summary, for all values of k and t the boundary conditions are either of limit-circle 
or of limit-point type. Furthermore, both end points are non-oscillatory and independent of 
the energy eigenvalue E. This allows us to apply theorem 7.5 in [30] which ensures that the 
spectrum is purely discrete, bounded below and unbounded above. 

If limit-circle endpoints are present, a Sturm-Liouville problem is not self-adjoint unless 
boundary conditions are imposed. For a non-oscillatory endpoint a there is always one solu- 
tion which is "small", known as the subdominant solution. More precisely, the subdominant 
solution w(A), unique up to a scalar factor, satisfies 

u(X) 

lim -U- = 0, (5.23) 

A->a v(X) 

where v (A) is any linearly independent solution, see theorem 7.15(i) in [30]. In our case, the 
subdominant solution is the non-singular solution, e.g. ciX 1 as A — > 0. A natural boundary 
condition is defined by the subdominance condition, i.e. always taking the non-singular solu- 
tion. Theorem 7.21 in [30] then guarantees that the subdominance condition defines a valid 
self-adjoint problem known as the Friedrichs extension. 
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5.4 Numerical results 



In the following, we briefly sketch our numerical scheme for Hq with k and t fixed. The 
equations for H\ and H 2 require minor, but straightforward modifications. We are using a 
multi-component shooting method which solves a collection of initial value problems (for the 
eigenvalue equation at fixed E) numerically with a standard adaptive Runge-Kutta method. 
The first initial value problem has initial values at Ao ~ and the solution a(A) of the 
differential equation (5.11) is evaluated at a matching point A m ^> A . Allowing the initial 
data to span the m = 2 min(/c, t) + l dimensional space specified by the boundary condition at 
0, this produces a 2m x m matrix whose columns consist of a(A m ) above a'(A m ). The second 
initial value problem has initial values at Aoo 3> A m and is again evaluated at A m , resulting, 
as the initial data span the boundary condition, in a 2m x m matrix 77 constructed similarly. 
By construction, a general solution satisfying the boundary condition at lies, at A m , in the 
column span of <ft, while a general solution satisfying the boundary condition at 00 lies, at A m 
in the column span of 77. Hence, E is an eigenvalue if these spans have nontrivial intersection, 
that is, if 

d(E) = det(0/7) = 0. (5.24) 

Having constructed d(E) numerically, we find its roots using a bisection method. Typical 
values of our constants are Ao = 0.001, A m = 3 and A^ = 30. 

Table 1 shows the lowest energy levels for Hq = |A which corresponds to the Laplacian 
for the L 2 metric, for Hi = |A + which also includes the curvature corrections, and 
for H2 = + + c io which includes curvature corrections and Casimir energy. Energy 
levels are rounded to two decimal places. The energy levels are ordered according to their 
respective energy, and are labelled by isospin quantum number k, total angular momentum 
quantum number t and parity P. Curly brackets indicate that the states (A;, t) and (t, k) 
have the same energy. States with the same quantum numbers but different energies form 
sequences of radially excited states. We have chosen to display levels with an energy not 
greater than the second excited (k — 0, t — 0) state, which are 19 energy levels in total. The 
columns "degeneracy" give the number of different states with the same energy. For k = t, 
the degeneracy is given by (2k + l) 2 , whereas for k 7^ t it is 2(2k + l)(2i + 1), where the extra 
factor of 2 arises from the (k, t) — > (t, k) symmetry of the spectrum. 

The spectra of Hq and H\ are remarkably similar. The order of the energy levels remains 
the same apart from two exceptions which are marked with a * in table 1 and which will be 
discussed in more detail later. The Casimir energy leads to significant changes, nevertheless 
the spectrum still shares some similarities. The ground state of H has energy ^ 0) = 0.00. 
The curvature term increases the energy of the ground state of Hi to E^ = 0.22, whereas the 
renormalised Casimir energy leads to a decrease in the ground state of H 2 to Eq = —0.41. 

In order to compare the energies of the excited states, we shift the spectrum of the ground 
states of Hi and H2 to 0.00 and then calculate the relative difference of the corresponding 
energy levels. Hence, the relative difference in energy between the nth excited state En^ of 
Hq and the nth excited state E^ of Hi is given by 



E^ - (e?> - E { 



(5.25) 
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Table 1: This table shows the lowest energy eigenvalues of Ho, Hi and if 2- Each energy 
eigenvalue is labelled with its degeneracy and its quantum numbers k and t and parity P. 
Energy eigenvalues marked with a * occur in a different order than their counterpart in the 
H spectrum. 



The relative difference of the {0, 1} _ states of Ho and H\ is 15%, and 7.7% for the first excited 
{0, 0}+ state. All the remaining states have a relative difference of less than 7%. 

As mentioned earlier, some energy levels change order. Therefore, it is useful to calculate 
the relative error of states with the same quantum numbers {j, t} p . The first transposition 
occurs for the first {0, 2}+ state and the first {1, 1}~ state. The relative difference between the 
energy of the {0, 2}+ state between H and Hi is 10.4%. The second transposition occurs for 
the first {0, 3}~ and the first even {1, 2}+ state, and the {0, 3}~ state has a relative differences 
of 8.0%. This leads to the observation that the difference between positive and negative parity 
states with the same {k, t} is reduced for Hi compared to H . Furthermore, Hi seem to favour 
states with k ~ t. 

As can be seen in table 1 the spectrum of H2 is shows many transpositions compared to 
the spectra of Hq and Hi. However, these transpositions occur for states which are close in 
energy, and the relative positions of the {0, 0}+ state and the {0, 0}+ excited states remain 
almost unchanged. Calculating the relative differences as in formula (5.25) shows that all 
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relative differences are less than 24%. 



5.5 Changing the radii 

It is interesting to consider how our results change if the radii R\, R 2 of the domain and target 
spheres are altered. We adopt the convention that a tilde signifies the case of general Ri, R2 
while undecorated variables refer to the case R\ — R 2 — 1. It is immediate from (2.4) that 
the L 2 metric on Rat n scales as 7 = R^R 2 /), and hence, from (3.2), we see that the Laplacian 
scales as A = (RiR 2 )~ 2 A. Now scalar curvature scales in the same way as the Laplacian 
(as can easily be seen, in this case, from the formula for k in [37], for example). Hence, the 
Hamiltonians Hq and Hi, scale homogeneously, Hi = (RiR 2 )~ 2 Hi, and so their spectra can be 
obtained from table 1 by a simple rescaling. Now the energies shown in this table are quantum 
corrections to the classical energy, which, by the Lichnerowicz bound is 4irR 2 . Hence, for Hq, 
Hi, the total energy of the kth energy level is 

Ef tal = ^Rl + -§^ 2 (5.26) 

where E^ is the eigenvalue of Hi. It is interesting to note that the quantum correction becomes 
(naively) dominant in the case of small target space (R 2 small). 

The behaviour of H 2 is more subtle because the Casmir energy scales differently from 
the other terms. To see this, we determine how the Jacobi operator for a harmonic map 
(j) : M —> N scales under homotheties of M and N. 

Proposition 8 Let </? : (M m ,g) — >■ (N n , h) be harmonic with Jacobi operator J? . Ifg = R\g, 
h = R\h, where Ri,R 2 > are constants, then the Jacobi operator of tp as a harmonic map 
(M,g) (N,h) is _ 

Proof: Let ip s> t '■ M — > N be a smooth two-parameter variation of if = (fo,o, with d s ip Si t \ s =t=o = 
X, d t (p St t\s=t=o = Ye Y^ip^TN). We have, in obvious notation, 

ym-2 D 2 



E(ip ait ) = R?- 2 R 2 2 E(<p St 
fc(X,F) = ^£ (¥>.,, 



= R?- 2 R 2 2 Hess(X, Y) 

s=t=0 



R[ l - 2 R\ [ h(X,/Y)vo\= [ h(X,Rf JY)vo\. 

J M JM 



□ 



Since ^ is (formally) a sum of square roots of eigenvalues of J? , we see that it scales as 
^ = R~[ Xc € . Hence, the scaling behaviour of H 2 is inhomogeneous, and, except in the special 
case that R\R 2 2 = 1, the spectrum of H 2 cannot be deduced directly from the spectrum of H 2 . 
Since ^ is independent of the radius R 2 whereas the spectra of H and Hi scale like 1/R 2 , the 
effect of the Casimir energy becomes less important for small R 2 . We have also numerically 
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checked that a smaller radius of the target sphere reduces the effect of the Casimir energy term 
significantly. For example for Ri = 1 and R2 — |, there is only one transposition, compared 
to the Hi spectrum. 

Throughout this paper we have used the convention that h = 1. If we reinstate the param- 
eter h then the classical energy scales like h°, the Casimir energy scales like h 1 whereas the 
Laplacian and curvature terms scale like h 2 . Hence, h can be removed by redefining P x y Ri/h, 
so the semi-classical limit K — > coincides (formally) with the planar limit Ri — > 00. 




0.2 - 
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s 

Figure 3: Energy levels interpolating between the L 2 metric and the CP 3 metric. Here solid 
lines correspond to positive and dashed lines to negative parity. 



5.6 Deformation to the Fubini-Study metric 

As a nontrivial test of our calculations we calculate how the spectrum changes as the metric is 
smoothly deformed from the L 2 metric to the well-known Fubini-Study metric of CP 3 , where 
the spectrum and degeneracy of the Laplace operator have been calculated explicitly. We 
consider the one-parameter family of SO (3) x SO(3) invariant kahler metrics defined, as in 
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Proposition 1, by the functions 



A S (X) = 32sA FS (\) + (1 - s)A L 2(X), 



< s < 1, 



(5.27) 



where Aps and A L i are given in (2.14) and (2.12), respectively, and the factor of 32 ensures 
that the eigenvalues of the two Laplacians are of the same order of magnitude. The pth 
eigenvalue of the Laplacian for the Fubini-Study metric on CP 3 with constant holomorphic 
sectional curvature 4 (and hence coefficient function Aps) is 



for p > 1 [6]. Although it is not made clear in [6], the degeneracy of E = must be 1 by the 
Hodge isomorphism theorem, since CP 3 is connected, hence dimif°(CP 3 ) = 1 which equals 
the dimension of the space of harmonic 0-forms. Note that deg(P p ) for p > 1 is an integer 
divisible by 9. In figure 3 we show how the energy levels change as the parameter s in (5.27) 
is increased from (the L 2 metric) to 1 (the Fubini-Study metric). We distinguish between 
states with even and odd parity. The ground state is E = for both metrics. In figure 3 
we follow all the energy levels of table 1, and it can be seen how the different energy levels 
become degenerate at s — 1, the Fubini-Study limit. Our numerically computed energy levels 
agree with the exact result to at least two decimal places, when we take into acount the factor 
\ in Hq and normalize the Fubini-Study energy by a factor of 32, which arises from the factor 
32 in (5.27). 

Our numerical scheme at s = 1 does not find all eigenfunctions of A on CP 3 , because, 
except for the ground state, we impose a boundary condition which forces the wavefunction 
to vanish on c^Rati, the boundary of Rati at infinity. It is known [37] that c^Rati coincides 
with the image of the inclusion CP 1 x CP 1 CP 3 , ([ao, ai], [bo, bi\) h-> [ao&i, ai&i, ai&i, ai^o]- 
The stabilizer of this subset of CP 3 in S77(4), the isometry group of CP 3 , is isomorphic to 
SU(2) x SU(2) (the isometry group of Rati itself). So, given an eigenfunction vanishing 
on dooRatx, the SU(4:) action generates a 9-dimensional orbit of degenerate eigenfunctions 
which do not vanish on dooRati (because dim SU (4) = 15 and dim SU (2) x SU(2) = 6). 
Hence, we expect the degeneracies found by our numerics (at s = 1) to be a factor 9 smaller 
than the degeneracies of the true Laplacian on CP 3 . This is precisely what we find in our 
numerics. For example, the degeneracy of the second CP 3 eigenvalue is deg(Pi)/9 = 15, which 
corresponds to the levels {0, 1} and {1, 1} with degeneracy 2*3 + 3*3 = 15, see table 1. 
Similarly, deg(P 2 )/9 = 84, corresponds to {0,0}, {1, 1}, {0,2}, {1, 1}, {1,2}, and {2,2} with 
1 + 3*3 + 2*5 + 3*3 + 2*3*5 + 5*5 = 84. Further, deg(P 3 )/9 = 300, which corresponds to 
{0, 1}, {1, 1}, {1, 2}, {0, 3}, {1, 2}, {2, 2}, {1, 3}, {2, 2}, {2, 3}, and {3, 3}, whose degenercies 
add up to 300, as expected. To display all the levels which contribute to the fourth energy 
eigenvalue of CP 3 we have also included the {3, 3} level in figure 3. Note that the second 
excited {0, 0} state contributes to the fifth energy eigenvalue of CP 3 . 



E p = 4p(3 + p) 



(5.28) 



with degeneracy 



deg(P p ) = f(3 + 2p)(p+l) 2 (p + 2) 2 



(5.29) 
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6 Concluding remarks 



In this paper we discussed the semi-classical quantization of soliton dynamics for CP 1 lumps 
moving on a 2-sphere. We followed Moss and Shiiki [28] who derived a Born-Oppenheimer 
approximation to the quantum dynamics based on the moduli space approximation. We 
were able to evaluate three different truncations of the Born-Oppenheimer Hamiltonian P#o? 
namely the geometric Laplacian H = ^A, the first geometric correction Hi given by the 
scalar curvature k of the moduli space, and the Hamiltonian H2 which consists of Hi together 
with the Casimir energy. The Casimir energy is notoriously difficult to evaluate. At A = and 
A = 00 the spectrum of the Jacobi operator is known explicitly. Our approach is to calculate 
the regularized Casimir energy using zeta- function regularization for A = and A = 00. Then 
the intermediate values are calculated using an approximate self-similarity of the eigenvalue 
spectrum. We calculated the first 19 energy levels for Hq, Hi and H2, and found that the 
first two spectra are remarkably similar whereas the inclusion of the Casimir energy leads to 
significant changes. There is an overall shift of the spectrum of Hi and H2. The relative errors 
between H and Hi are 15% for the first excited state and less that 8% for other excited states. 
By contrast, the relative errors between H Q and H 2 are as high as 24% for some excited states. 
There are only two transpositions of states where energy levels of Hi do not have the same 
order compared to the H spectrum. However, the spectrum of H 2 shows many transpositions. 

We proved that all the spectra enjoy a rather sursprising spin-isospin exchange symmetry. 
The proof rested on the identification of a hidden isometry of Rati which becomes manifest only 
after one identifies Ratx = TSO(3). As a non-trivial check of our calculations we calculated 
how the energy levels change when we interpolate between the L 2 metric and the Fubini-Study 
metric on CP 3 . We reproduced the known spectrum for CP 3 and discussed the degeneracies 
of the energy levels. 

Our approach allows us to calculate the spectrum of the Laplacian for any 5*0(3) x SO (3) 
invariant kahler metric on Rati. One interesting choice is A = c/A, for which the coefficient 
of (A • T) 2 vanishes in (5.1). In this case, the angular momentum operator L 2 also commutes 
with the Hamiltonian and I becomes an additional quantum number. 

The lump dynamics serves as a toy model for other solitons. For example in three spatial 
dimensions, two physically relevant models are the Skyrme model [34] and the Faddeev-Hopf 
model [10, 11]. In both models, the solitons can be quantized as fermions due to so-called 
Finkelstein-Rubinstein constraints [12]. Adkins, Nappi and Witten first quantized the B = 1 
Skyrmion in [2]. The effects of the Casimir energy on the predictions of the Skyrme model in 
the B = 1 sector have been discussed in great detail in [26]. The authors calculated the 1-loop 
corrections to various physical quantities using phase-shift techniques to evaluate the Casimir 
energy. For higher topological charge the ground and lowest excited states in the Skyrme 
model and the Faddeev-Hopf model were calculated in [15, 18, 4] and [19] using zero-mode 
quantization and Finkelstein-Rubinstein constraints. The results of this paper suggest that the 
order of the states would remain the same, with minor changes, if the curvature correction in 
the Born-Oppenheimer approximation was taken into account. However, the Casimir energy 
could lead to significant changes in the order of states. A more careful analysis of higher order 
terms in these models, and in particular, of the Casimir energy, would be very useful. 
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Appendix 

A The rough Laplacian 

In this appendix we will compute in detail the action of the rough Laplacian on sections 
Y of p~ l TN of the form 

Y = a(9)cosm<j)E 1 . (A.l) 

Recall that (M, g) = (N, h) = the unit sphere, (9, (ft) are the usual polar coordinates on M or 
N, Ei = d/d8, E 2 = cosec8d/d<ft is an orthonormal frame on M or N, E\ = E\ o tp, E 2 o <p, 
and tp : M — > N is a map of the form tp(9, (ft) = (f(9), (ft) in polar coordinates. Now, given a 
map between Riemannian manifolds p> : M — > N, the pullback connexion is the unique 
connexion on the vector bundle p~ l TN satisfying the axioms 

V x {Y l + Y 2 ) = V£yi + V£y 2 (A.2) 
V x {fY x ) = Xim + fV^ (A.3) 
V x (Yop) = (V^ x Y)op (A.4) 

for all X E r(TM), Y U Y 2 E T(p~ 1 TN), f E C°°(M) and Y E T(TN), where is the 
Levi-Civita connexion on (N, h). In our case (which coincides with V M ) is determined by 
its action on the frame {£1, E 2 }, 

V N E 1 = co%9e 2 ®E 2 , V N E 2 = -cot0e 2 ® E u (A.5) 

where {e\,e 2 } is the coframe dual to {Ei,E 2 }. From this we may deduce how V 1 ^ acts on 
Et,E 2 . By property (A.4) and (A.5), 

(A.6) 
(A.7) 

(A.8) 
(A.9) 

The rough Laplacian has 4 terms, 

= -V^(V|F) - V| 2 (V| Y) + % ME Y + V% UE Y. (A.10) 
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^e 2 E 2 = 


sin 6 a2 


O (p = 


cos f ~ 

-— ^1 
sm 9 
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We evaluate these in turn, for Y = a(6) cosmipEi. First 

d6 2 



d 2 ~ 

V^V^acosm^Ei) = ——(acosm(f))Ei = a" cos m<f>Ei. (A. 11) 



The second term is more involved: 



V^ 2 (acosm0i?i) = aE 2 [cos m<p] E\ + a cos m0V^ 2 E\ 

sin md) ~ cos md> cos f ~ , . _ . 

= -ma— -p£i + a £ 2 (A.12) 

sin t7 sin # 

^ , ^ ~, 2 cosm0~ sinm^cos/- 

sin 6 1 sin 6 1 



cosm0cos 2 / ~ 

" a ^7; ^i- (A. 13) 

sin & 



The third term vanishes since V^-Ei = 0, while the fourth term is 

V^ M „ (acosm<j)Ei) = — cot^V^ (acosm0i?i) = — a cot # cosm<f>Ei. (A. 14) 

Assembling the pieces, one sees that 



~ / 7T? COS 2 f \ ~ 

Au,(a cos m<f)Ei) = I —a" — a' cot 6> H 5— a H 5 — a I cos m<pEi 

\ sin sin - 9 J 

( Ira cos / \ ~ 
+ — , n - a sinm0£/ 2 , (A. 15) 

\ sm / 

as claimed in equation (4.14). 

B Angular momentum calculations 

In this appendix we briefly describe the evaluation of the operator (A • T) 2 . We closely follow 
the notation in [8]. Note the spin quantum number j satisfies j = k since J 2 = K 2 . We choose 
the following convention for the spherical harmonics, 



(2/+l)(/ m)! 
47r [i + my. 



where the associated Legendre polynomials are given by 



P™{x) = 1 X J g— (x 2 - I) 1 . (B.17) 



Note that Y* = (-l) m Y, _ m and 



Urn = r l,-m 

2tT TV 
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It is easy to show that (A • T) = (A • J), and we showed in section 3 that this operator 
commutes with T and J 2 . Hence, the operator does not change the quantum numbers t, k, k% 
and £3. Furthermore, the matrix elements are independent of £3 and k 3 , hence we can choose 
without loss of generality £3 = and k% = 0. We need to evaluate 

M = (M h ) = (k, l,t,0\\ ■ J\k,l, t, 0), (B.19) 

which is a (2 mm(k, t) + l) by (2 mm(k, t) + l) matrix, because the angular momentum quantum 
numbers satisfy min(/c, t) < I, I < k + 1. Since A • J is a Hermitian operator the matrix M is 
Hermitian, namely 

M U = M\. (B.20) 
Using Clebsch-Gordan coefficients (B.19) can be rewritten as 

M = (k, I, t,0\k, mi, I, —mi)(k, mi, I, — mi\ A ■ J|fc, m2, 1, —m2){k, m2, 1, —m2\k, l,t,0). 



mi ,m,2 



(B.21) 



The operator A ■ J can be written in terms of spherical harmonics as 



Hence, we only need to evaluate 

(k,m 1 \J + \k,m 2 ) = S mi>m2+1 ^ (k - m 2 )(k + m 2 + 1), (B.23) 

(k,mi\J_\k,m 2 ) = <W,m2-i\/P + m 2 )(k - m 2 + 1), (B.24) 

(k,mi\J 3 \k,m 2 ) = 5 mum2 m 2 , (B.25) 

and 

(l,mi\Y l:m \l,m 2 ), (B.26) 
which can be calculated using (4.6.3) in [8], namely, 

2-7T 7T 

Y lumi (6, <f))Y l2 , m2 (6, <f>)Y l3:Tn3 (6, <P) sin 9d0d<f> (B.27) 








(2/ 1 + l)(2/ 2 + l)(2/ 3 + l) (li l 2 l 3 \ / Zi Z 2 h 
A-k \ J \ mi m 2 m 3 

where the two matrices in (B.28) correspond to 3j symbols. The 3j symbol 

I I 1 




(B.28) 



(B.29) 



vanishes unless Z = Z ± 1, which can be used to show that M vanishes in the cases k = or 
t — 0. For example, for k = the only allowed values of the angular momentum are I = t and 
I = t, which vanishes since I 7^ I ± 1. 
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Following the notation in [8], the Clebsch-Gordan coefficients and the 3j symbols have 
real entries. Hence, M is a real matrix and (B.20) implies that M is symmetric. Note 
that the eigenvalues of M are 0, ±1, ±2, . . . , ±min(/c,t), which can be explained as follows. 
Heuristically, (A ■ T) is the projection of the total angular momentum operator onto the A 
direction. Hence, it can be rotated by a change of variables to T3 which has eigenvalues 
0, ±1, • • ■ ± min(/c, t). The highest possible value of £3 is rnin(fc, t) because the dimension of M 
is 2 min(fc, t) + 1. 

Now, we can evaluate the matrix (B.19) for various values of k and t using Maple, for 
example. For small values the formulae are more tractable. For example for (k,t) = we 
obtain 



We showed in section 5.1 that there is a (k,t) —> (t,k) for the respective matrices M. This 
was verified by explicit calculation for all values (k,t) in table 1. As can be seen in (B.30) 
there is a natural chess board structure in the matrix M 2 , which follows from the fact that M 
vanishes unless I = I ± 1, and hence M 2 vanishes unless / = Z + 2, Z, or/ — 2. 
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